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The "FT-symmetric Hamiltonian H — p 2 + x 2 (ix) £ (e real) exhibits a phase transition at e — 0. 
When e > 0, the eigenvalues are all real, positive, discrete, and grow as £ increases. However, when 
£ < there are only a finite number of real eigenvalues. As £ approaches —1 from above, the number 
of real eigenvalues decreases to one, and this eigenvalue becomes infinite at £ = — 1. In this paper 
it is shown that these qualitative spectral behaviors are generic and that they are exhibited by the 
eigenvalues of the general class of Hamiltonians H^ 2n ^ — p 2n -\-x 2 (ix) £ (£ real, n = 1, 2, 3, . . .). The 
complex classical behaviors of these Hamiltonians are also examined. 
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I. EIGENVALUE BEHAVIORS 

The well-studied family of PT-symmetric Hamiltoni- 
ans 

H = p 2 + x 2 (ix) £ (e real) (1) 

is a continuous complex deformation in the parame- 
ter e of the conventional harmonic-oscillator Hamilto- 
nian p 2 + x 2 . At e = the eigenvalues of H are 
the usual harmonic-oscillator eigenvalues E n = 2n + 1 
(n = 0, 1, 2, ...), but as e varies away from 0, the 
eigenvalues exhibit several characteristic qualitative be- 
haviors: When £ is positive, the eigenvalues are all real, 
positive, and discrete, and the eigenvalues increase as e 
increases. When e is negative, most of the eigenvalues 
are complex. As e decreases from to —1, the eigenval- 
ues sequentially become pairwise degenerate and enter 
the complex plane as complex-conjugate pairs. As e con- 
tinues to decrease, eventually only one real eigenvalue 
remains, and this eigenvalue becomes infinite at e = — 1 
[6]. This behavior is displayed in Fig. [I] 

The parametric region e > is called the region of un- 
broken VT symmetry because all of the eigenfunctions of 
the Hamiltonian are simultaneously eigenstates of VT] 
the parametric region £ < is called the region of broken 
VT symmetry because the eigenfunctions of the Hamilto- 
nian are not all eigenstates of VT [7] . A phase transition 
at £ = separates the regions of unbroken and broken 
symmetry. Parametric families of T^T-symmetric Hamil- 
tonians typically exhibit a phase transition with entirely 
real eigenvalues on one side and complex eigenvalues on 
the other. The VT phase transition has been observed 
repeatedly in laboratory experiments [5H22]. 
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FIG. 1: Energy levels of the parametric family of Hamiltoni- 
ans |l]). When £ > 0, the eigenvalues are all real and positive, 
and increase with increasing e. When £ decreases below 0, 
the eigenvalues disappear into the complex plane as complex 
conjugate pairs. Eventually, only one real eigenvalue remains 
when £ is less than about —0.57, and as £ approaches —1 from 
above, this eigenvalue becomes infinite. 

The purpose of this paper is to report our observation 
that the qualitative features of the eigenvalues of H in 
([!]), as shown in Fig. [I] are generic; that is, the same 
qualitative behavior is displayed by the general family of 
Hamiltonians 

H {2n) = p 2n + rf^y ^ ^ n = i? 2 , 3, . . .) (2) 

for which H in ([!]) is just . For example, when n = 2, 
the time-independent eigenvalue problem H^ip = Eip 
in coordinate space becomes the fourth- order differential 
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FIG. 2: Eigenvalues of the T^T-symmetric Hamiltonian 

in |2| plotted as functions of e for — 1 < e < 4. The eigenvalue 

behavior agrees qualitatively with that in Fig. [I] for the case 



FIG. 3: Eigenvalues of the VT- symmetric Hamiltonian 
in j2]) plotted as functions of e for — 1 < e < 2. The eigenvalue 
behavior is qualitatively similar to that in Figs, [I] and [2] for 
tf< 2 > and i/W. 



equation 



dx 4 



ip{x) + x 2 (ix) e ip(x) = Eip(x). 



(3) 



To determine the boundary conditions on ^(x), we begin 
with the requirement that ip(x) vanish exponentially fast 
as \x\ — >• oo on the real-x axis when e = 0, so the eigenval- 
ues are the same as those of the conventionally Hermitian 
quartic-anharmonic-oscillator Hamiltonian p 2 + x 4 [23] . 
Then, this eigenvalue problem is extended into the com- 
plex domain by allowing the parameter e to range con- 
tinuously from —1 to +oo. As we explain in Sec.[Hj this 
analysis determines the boundary conditions on ip{x) in 
the complex- x plane for real — 1 < £ < oo. 

The eigenvalues of are plotted in Fig. [5J The 

behaviors of the eigenvalues in Fig. [2] are qualitatively 
similar to the eigenvalues in Fig. [T] The eigenvalues in 
both cases are real, positive, and discrete for e > and 
they become complex below the phase transition at e = 0; 
at e = — 1 there are no eigenvalues at all. 

When n = 3, the time-independent eigenvalue problem 
H^ip = Eip is the sixth- order differential equation 



dx 6 



ijj(x) + x 2 {ix) £ ip{x) = EtJ)(x). 



(4) 



When e = 0, the eigenvalues are the same as 
those of the conventionally Hermitian sextic-anharmonic- 
oscillator Hamiltonian p 2 + x 6 . Observe that for —1 < 
£ < 00 the eigenvalues, which are shown in Fig.[3J behave 
qualitatively like those in Figs, [l] and [2] 



This paper is organized very simply: In Sec.|ll|we ex- 
plain the asymptotic and numerical analysis that leads 
to the results presented in Figs. [2] and [3] In Sec. HI 



study the complex classical trajectories associated with 
the Hamiltonians H^ 2 \ H^ A \ and for a variety of 
values of e. These trajectories demonstrate some qual- 
itative differences between these Hamiltonians. Finally, 



in Sec. IV we make some brief concluding remarks. 



II. ASYMPTOTIC AND NUMERICAL 
ANALYSIS 

In this section we explain the techniques for calculating 
the eigenvalues that are graphed in Figs. [T}j3] We begin 
by reviewing the procedure for finding the eigenvalues of 

//•-•. 



A. Review of the case 

The properties of the eigenvalues of were first 

elucidated in Ref. [1 . For completeness we give a brief 
review of the techniques that were used to construct an 
accurate plot of the eigenvalues. 

The eigenvalue equation (the time-independent 
Schrodinger equation) for reads 



il)"{x) + x 2 (ixf^(x) = Eil){x). 



(5) 



This differential equation has many different sets of eigen- 
values, which are distinguished by the boundary condi- 
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tions imposed on the eigenfunctions. Thus, it is crucial 
to identify the specific boundary conditions on the eigen- 
functions ip(x) that give rise to the eigenvalues plotted 
in Fig. [I] 

WKB theory can be used to obtain the possible asymp- 
totic behaviors of ip(x) for large \x\: 

i/j(x) ~ exp (±x 2+£ /V" £ / 4 ) (\x\ -> oo), (6) 

where we have used i = e 27r / 2 . Since ([5]) is a second-order 
differential equation, there are two solutions; we want 
the solution that decays exponentially as \x\ — >• oo on the 
real-x axis when s = 0. 

To understand what happens when e ^ 0, we let x = 
re %d and require that as r — » oo, 



at r = ro, we use the WKB approximation for large tq\ 



ip(r) ~ exp 



_ r 2+ £ /2 e i[0(2+ £ /2)+7r £ /4] 



Note that ip(r) vanishes in a Stokes wedge. At the center 
of the this wedge the phase, which is the expression in 
square brackets in (|7|), vanishes. This gives the angular 
location of the center of the wedge: 



9 R = -7rs/(8 + 2e). 



(8) 



We use the notation Or here to indicate that this is the 
right Stokes wedge. There is also a left Stokes wedge, 
which is the mirror image of the right Stokes wedge; the 
angle of the center of the left Stokes wedge is called 0l- 
The mirror for this problem is the imaginary axis; the 
imaginary axis serves as the mirror because the Hamilto- 
nian is VT symmetric. Note that Or = at e = 
and that Or — » — tt/2 as e —> oo. 

The opening angle AO of the right Stokes wedge (and 
of the left Stokes wedge) is determined by requiring that 
the phase in square brackets in ^ is ±7r/2. This gives 



AO = 2tt/(4 + £). 



(9) 



The opening angles of the Stokes wedges shrink to as 
e — >• oo. However, as e —> — 1, the opening angles expand 
to AO/2 = 7r/3. Thus, since Or = tt/6 at e = — 1, we 
see that the upper edge of the right wedge is at Or + 
AO/ 2 — 7r/2. At this value of e the left and right wedges 
meet and the eigenvalue problem disappears because the 
integration contour can be deformed to infinity. This 
is an intuitive way to understand the rigorous result of 
Herbst that the spectrum becomes null when e = — 1 [6]. 

The numerical calculation of the eigenvalues uses the 
standard shooting method. We integrate numerically with 
respect to the r variable down the centers of the left and 
right wedges and join the two solutions at the origin. The 
differential equation in the right wedge in the r variable 
reads 

^k(r) = {r 2+£ + Ee™«) ^(r) = Q(r)V*(r). (10) 

The numerical integration begins at a large value of r, say 
r = ro- To determine the appropriate initial conditions 



^wkb Oo) ~ KQ 1/4 (r ) exp 



dry/Q(r) 



(ii) 



where K is an arbitrary constant. We then adjust the 
constant K so that 

ipR M = 1 and i>' R (r ) = - y/Q{r ). (12) 

For the eigenvalues shown in Fig. [l] we choose ro = 12. 

We patch the left and right solutions at the origin as 
follows: In the left wedge the eigenfunction that vanishes 
exponentially for large argument s is 0^(5), and in the 
right wedge the eigenfunction that vanishes exponentially 
for large r is f3ipR(r), where a and (3 are arbitrary con- 
stants. Thus, the patching conditions are 

a^L(O) = £<M0), ae-^^(O) = ^"^^(O)- (13) 

These are homogeneous linear equations, and thus by 
Cramer's rule the determinant D must vanish, where 



D = det 



-**^(o) 



-Mo) 
~ iGL Mo) 



(14) 



The shooting procedure consists of adjusting the eigen- 
value E so as to minimize \D\. This numerical procedure 
is highly efficient and gives a rapidly convergent sequence 
of numerical approximations to E. 



B. Case if (4) 

Let us consider the fourth-order differential-equation 
eigenvalue problem in (|3|. A local analysis of this differ- 
ential equation for large \x\ shows that, apart from mul- 
tiplicative algebraic corrections, the four possible asymp- 
totic behaviors of ip(x) are 



ip(x) ~ exp 



4uj 
6 + 6 



, x (6+£)/4 e i7r £ /8 



(M 



(15) 



where uo A = —1. We must exclude two of these four 
solutions because they grow exponentially when \x\ ^> 1. 

We now describe the Stokes wedges that are associated 
with the solutions. We let x = re z0 , and in the right 
Stokes wedge we set = Or. Neglecting for the moment 
the effect of cj, the phase angle in (15) is (6 + £)0r/A + 
tts/8. This angle vanishes at the center of the Stokes 
wedge. Thus, 



Or 



-7re/(12 + 2e). 



(16) 



As in the case of H ^ , the center of the wedge is located 
at Or = when e = 0, and the center rotates downward 
(clockwise) to — tt/2 as e — >• 00. 

As before, the opening angle A# of the Stokes wedge is 
determined by the condition that the phase has the value 
±7r/2. This gives 



A(9 = 4tt/(6 + £:). 



(17) 
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Thus, the opening angle is 2ir/3 at e = and the opening 
angle vanishes as e — » oo. 

We now observe one of the universal features of the 
Hamiltonians in At £ = -1, Or = n/lO and A (9/2 = 
47r/10. Thus, the upper edge of the right wedge is at 
R + A0/2 = tt/2. By symmetry, there is a left 
wedge (a mirror reflection of the right wedge) in the left- 
half plane. The center of the left wedge is at the angle 
Ol = — 7r + Or. Hence, when e = — 1, the left and the 
right wedges touch, the eigenvalue problem collapses, and 
the spectrum becomes null. (This is the generalization 
of the effect that Herbst found for the Hamiltonian H = 
p 2 ±ix [6].) This collapse of the eigenvalue problem is a 
universal phenomenon and always occurs at e = —1 for 
all Hamiltonians in ([2|. 

We now explain the numerical procedure for finding the 
eigenvalues in Fig. [2] In the right wedge we substitute 
x = re %dn in the differential equation (p3|) and obtain 



v4 4) (r) = -e 



Ee 



(18) 

There are two solutions that grow exponentially with r 
and two solutions that vanish exponentially with r. 

Until now, u has not played any part in the analysis. 
However, it plays a crucial role in determining the bound- 
ary conditions. We integrate numerically downward in r 



from a large value of r = r$. To fix the initial condi- 
tions we use the WKB approximation for fourth-order 
differential equations [24 : 



^wkb (r ) ~ KQ s ( ro ) ex p 



pro 

/ drQ^ 4 (r) 
Jo 



(19) 



where u 4 = —1. We choose the constant K so that 

^(r ) = l- (20) 

Then, using uj\ — e 37 ™/ 4 and c^2 = e -37 ™/ 4 so that the 
solution decays exponentially, we take 



Vffi-M = ^Q W4J (r ) 



= ,.,*n(*/ 4 ) 



(21) 



(k) 

where ip y R j is the fcth derivative with respect to r. We 
then choose a large number for ro, say ro = 12. The 
solution in the left wedge is 



and the solution in the right wedge is 



(22) 



(23) 



We must now patch the left and right solutions at r = 0. To do so we require that the eigenfunctions and ^r 
be continuous and have continuous first, second, and third derivatives at r = 0. These four conditions are written as 
a system of four linear equations for a±, oli, and fc: 



«l^L,l(0)+a 2 ^L,2(0) = 01^fl,l(O)+ 02^,2(0), 



OL\e 
ot\e 
ot\e 



2ie ^L,i(°) + a 2 e- 2i ^L, 2 (0) = /3ie- 2 ^^ a (0) + /3 2 e- 2i ^ 2 (0), 
m ^lM+oc2e- m ^lM = /3ie- 3 ^< a (0)+^e- 3 ^< i2 (0). 



(24) 



Because these patching conditions are homogeneous equations, they translate via Cramer's rule to the condition that 
the following 4x4 determinant of the coefficients vanishes: 



D = det 



/ ^,i(0) 
e-**Vfl,i(0) 
-- 2iflB VSi(0) 



^, 2 (0) 
-^^,2(0) 

-2iS R 



„-3ifl 



l ^',i(0) 



-3i8 R 



^, 2 (0) 
^ 2 (0) 



-^L,l(0) 

-e-^^,i(°) 
-e- 2i ^Li(0) 



-^L,2(0) \ 

-e-^^, 2 (0) 



— e 
— e 



- 2 ^^ )2 (o) 

-3i6 Lq pl/l 



L,2 



(25) 



(0)/ 



The objective of the numerical shooting procedure is to 
find the values of E that minimize this determinant. 



tions) are given by 



600 
8+V 



lP(x) ~ exp ( _^ x (8+s)/6 e i7rs/12 



>h4 ( 2 6) 



C. Case if (6) 

For the eigenvalue problem in Q the possible asymp- 
totic behaviors for large \x\ (excluding algebraic correc- 



where uj = 1. 

There are six solutions, but we seek the three solutions 
that vanish as \x\ —> 



^ oo. We substitute x = re z ° R in 
(26), and we require that the solution vanish as r — » oo. 



Apart from cj, the phase angle in the asymptotic behavior 
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vanishes when 



9r = -7re/(16 + 2e). 



(27) 



This determines the center of the right wedge. At e = 0, 
0r = and as £ — >• oo, 0r — >• — 7r/2, just as before. 

The opening angle A# of the wedge is determined by 
requiring that the phase have the value ±tt/2. This gives 



AO = 6tt/(8 + £). 



(28) 



At e — the opening angle is 37r/4 and this angle vanishes 
as e —> oo. 

From ([27]) and ([28]) we see that at e = -1, R = 7r/14 
and A#/2 = 67r/14. Thus, the upper edge of the wedge is 
at 6 R + A0/2 = 7r/2. At this value of e, the left and right 
wedges touch and the eigenvalue problem disappears, just 
as in the case of and 

To find the eigenvalues using the numerical shooting 
method, we integrate down the center of the left and 
right wedges and patch the solutions at the origin. In 
the right wedge the differential equation in the r variable 



V4 6) (r) = (r 2+£ + Ee 6ie «) ^ R {r) = Q(r)^ R (r). (29) 

In each wedge there are three exponentially decaying so- 
lutions and three exponentially growing solutions. We 



select the decaying solutions by applying the appropri- 
ate boundary conditions at a large value of r = ro = 12. 
To identify the appropriate initial conditions for the nu- 
merical integration, we use the WKB approximation for 
large r : 



V'wkb (r ) ~ KQ 12 ( ro ) exp 



/ drQ^ 6 (r) 
Jo 



, (30) 



where ui = —1, co>2 = e 27 ™/ 3 , and ujs 
adjust the constant K such that 



Then, 



(ro) = 1. 



^](r )=^Q^ (r ). 



e- 2 ™/ 3 . We 
(31) 



The solution in the left wedge is 
and the solution in the right wedge is 



(32) 



(33) 



(34) 



The patching conditions give a vanishing determinant 
via Cramer's rule, where the determinant is 



(0) 



i(0) 



^R,l 

- 3 ^^', 1 (o) 
- 4ie ^S(o) 



-5i0 R „/,( 5 ) 



e-^^, 2 (0) 



V>fl, 2 (o) 

i0_ 

- 2i9 «< 2 (0) 

< 2 (0) 
: < } 2(0) 



o -3i0 
o —5i0R 



^,3(0) 
"^,3 



R,3 

-ue R „iSA) 



^ 5) 2 (o) 



-5i0R„/,( 5 



<Pr,3 



(0) 
(0) 
(0) 
(0) 
(0) 



— e 
—e 
—e 
—e 
—e 



-^L,l(0) 

-^^,i(0) 
2 ^^,i(0) 

3^^/(0) 
,(4) 



- 4 ^<i(o) 



—e 
—e 
—e 



— e 



— e 



-^L,2(0) 

-*^ i2 (o) 



'£, 2 (0) 
<£ 2 (o) 
^ 4) 2 (o) 

-5*e L ^,(5) 



-V>L, 3 (0) 

-e-^VL, 3 (0) 

-e-^V2. 3 (0) 
— e 

-4i<9 



(0) 



3 ^< 3 (o) 



— e 



—e 



^g(O) 

- Bi ^vS(o) y 



We shoot so as to minimize this determinant. 



III. CLASSICAL TRAJECTORIES 

While the eigenvalues of the quantum-mechanical 
Hamiltonians H^ 2n ^ are qualitatively similar, the 
classical-mechanical trajectories for these Hamiltonians 
exhibit marked differences. Using the numerical tech- 
niques described in Ref. [2] , we have studied the complex 
classical trajectories associated with the Hamiltonians in 
^ for e = 0, 2, 4, and 6. In Fig. [4] we plot four tra- 
jectories in x space and in p space for n = 1 for these 
values of e and in Fig. [5] we show corresponding plots for 
n = 2. Figures [6] and [7] display the corresponding plots 
for n = 3, and in these figures we can see the emergence 
of new kinds of topological structures. Figures [4][7| show 



a progressively increasing topological complexity, which 
arises because of the higher powers of p in the Hamiltoni- 
ans. It is interesting that the pattern of the eigenvalues 
in Figs. [l][3] shows no corresponding qualitative changes. 



IV. CONCLUDING REMARKS 

We have shown that the key qualitative features of 
the eigenvalue problem associated with #( 2n ) in ^ are 
universal. In particular, (i) the critical point between 
the regions of broken and unbroken symmetry is always 
at e = 0, (ii) the lower end of the broken region is always 
at e = —1 at which point the spectrum becomes null, and 
(iii) the eigenvalues are positive and rise with increasing 
s in the unbroken region. 

Qualitative differences between these Hamiltonians are 
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FIG. 4: Complex classical trajectories of the Hamiltonians 
if( 2 ) = p 2 -f x 2 (ix) e , where £ = (first row), e = 2 (sec- 
ond row), £ = 4 (third row), and e — 6 (fourth row). The 
left column is the projection of the classical trajectory in the 
complex- x plane and the right column is the projection of the 
trajectory in the complex-p plane. In every case the classical 
energy is 1. For the cases e — 2, 4, and 6 two trajectories are 
shown in the left column (blue and red online). The corre- 
sponding complex-p-plane trajectories do not cross, but rather 
lie on a two-sheeted Riemann surface. The p-plane trajecto- 
ries appear to be coincident because they are projected onto 
a single complex plane. 
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FIG. 5: Same as Fig.|4]panel except that n — 2. However, in 
this case we do not observe up-down trajectory pairing in the 
complex-x plane. 



much easier to spot at the classical level. We see that as 
n increases, the classical trajectories develop new and 
richer topologies. Near the turning points the trajecto- 
ries share the same features as those studied in Ref. [25] , 
in which the classical paths for higher-derivative Hamil- 
tonians were also investigated. 

CMB thanks the Leverhulme Foundation and the 
U.S. Department of Energy for financial support. Math- 
ematica was used to perform numerical calculations. 
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